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C^l ' The Generalized Relativistic Effective Core Potential (GRECP) method is described which allows 

one to simulate Breit interaction and finite nuclear models by an economic way and with high accu- 
racy. The corresponding GRECPs for the uranium, plutonium, eka-mercury (Eff2), eka-thallium 
fc^H ■ (E113) and eka-lead (Eff4) atoms are generated. The accuracy of these GRECPs and of the RECPs 

<H ' of other groups is estimated in atomic numerical SCF calculations with Coulomb two-electron inter- 

actions and point nucleus as compared to the corresponding all-electron Hartree-Fock-Dirac-Breit 
• calculations with the Fermi nuclear charge distribution. Different nuclear models and contributions 

of the Breit interaction between different shells are studied employing all-electron four-component 

1 i methods. 
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Investigation of physical and chemical properties of recently synthesized relatively long-living isotopes of superheavy 
elements (SHEs) with the nuclear charges Z=105 to 116 [1, 2, 3, 4] and their compounds is of fundamental importance 
. , for science. Their experimental lifetimes reach several hours now and the nuclei near the top of the "island of stability" 
are predicted to exist for many years. The experimental study of SHE properties is very difficult because of their 
extremely small quantities, only single atoms are available for research. Accurate calculations for SHEs and their 
O • compounds are needed in order to better understand their physical and chemical properties that often differ from 
'— —*[ those of the lighter homologs in the chemical groups due to very strong relativistic effects on their electronic shells. 
Besides, for elements decaying by spontaneous fission, the chemical identification is the only way to prove their Z 
number. 

j — . Experimental investigations of spectroscopic and other physical-chemical properties of actinides are severely ham- 
pered by their radioactive decay and radiation which lead to chemical modifications of the systems under study. The 
diversity of properties of lanthanide and actinide compounds is unique due to the multitude of their valency forms 
V} ' (which can vary over a wide range) and because of particular importance of relativistic effects. They are, therefore, 
of great interest both for fundamental research and for development of new technologies and materials. The most 
important practical problems involve storage and processing of radioactive waste and nuclear fuel, as well as pollution 
of the environment by radioactive waste, where most of the decayed elements are actinides. 

From the formal point of view, four-component correlation calculations [5, 6] based on Dirac-Coulomb-Breit (DCB) 
Hamiltonian (see [7, 8, 9, 10, 11] and references) can provide a very high accuracy of physical and chemical properties 
for molecules containing heavy atoms. However, such calculations were not widely used for such systems during last 
decade because of the following theoretical and technical complications [12]: 
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too many electrons are treated explicitly in heavy-atom systems and too large number of Gaussians is required 
for accurate description of the large number of oscillations, which valence spinors have in heavy atoms; 

the necessity to work with four-component Dirac spinors leads to serious complication of calculations as compared 
to the nonrclativistic case: 

(a) the number of kinetically-balanced two-component ("2c") uncontracted Gaussian basis spinors for the 
5mall components, Aj c , can be estimated as 2N^ C , where N^ c is the number of basis spinors for Large 
components; so the total number of uncontracted Gaussian basis spinors in the relativistic four-component 
("4c") calculations N^ s ~ 3A 2c and the number of two-electron integrals as [12] 



Ntmt ~ (l+2-2 2 +2 4 )iV 2 2 e c /nt = 25-N£ 
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Note, however, that the situation is seriously improved here during last years, see [5, 6, 13, 14]. 

(b) the number of basis 2c-spinors, Nj£ s , is twice more than the number of nonrelativistic basis one-component 
("lc") orbitals, N^ s , therefore 

The minimal number of two-electron integrals in the spin-orbit basis set, which are required to be saved 
coincides, obviously, with A^/nf 

The Relativistic Effective Core Potential (RECP) method is most widely used in calculations on molecules containing 
heavy atoms [15, 16] because it reduces drastically the computational cost at the integral generation, SCF and integral 
transformation stages. In our papers [17, 18, 19], the conventional radially-local (semi-local) form of the RECP 
operator (used by many groups up to now but suggested and first applied about 40 years ago [20, 21, 22]) was shown 
to be limited by accuracy and some nonlocal corrections to the RECP operator were suggested [17, 19, 23, 24], which 
have already allowed us to improve significantly the RECP accuracy [18, 19, 25, 26]. 

It is known that the Breit interaction can give contributions in excess of one thousand wave numbers even to 
energies of transitions between lowest-lying states of very heavy elements (see, e.g., tables IV and V). It is also clear 
that the point nuclear model becomes less appropriate when the nuclear charge is increased. Therefore, the RECPs 
designed for accurate calculations of actinidc and SHE compounds should allow one to take into account the Breit 
interaction and the finite size of nuclei. The most economic way is to incorporate the corresponding contributions 
into the RECP operator. 

I. GENERALIZED RECP METHOD 

In a series of papers (see [17, 18, 19, 24, 27] and references), we introduced and developed the Generalized RECP 
(GRECP) method. Its main features are: 

• The inner core (IC), outer core (OC) and valence (V) electrons are first treated employing different approxima- 
tions for each (including relaxation of the IC shells which are explicitly excluded from GRECP calculations). 

• GRECP involves both radially-local, separable and Huzinaga-type potentials as its components and particular 
cases. 

• The GRECP operator includes terms of other types for economical treatment of transition metals, lanthanidcs 
and actinides (see sections IA-IB). 

• The outer core pseudospinors (nodclcss) together with valence pseudospinors (nodal) are used for constructing 
the GRECP components [28]. 

• Quantum electrodynamics effects (see [29] and section IC), arbitrary nuclear models, and correlation with IC 
shells [30] can be efficiently treated within GRECPs. 

The GRECP method is described in detail in the above papers and we only add here that it allows one to avoid 
the complications of the four-component calculations described in the introduction (see also [30]) and to attain very 
high accuracy, limited in practice by possibilities of the correlation methods, while requiring moderate computational 
efforts when the IC, OC and V subspaces are appropriately chosen. 

The contributions of different nuclear models which are described by local potentials can be easily taken into account 
in the framework of the (G)RECP method. The situation is more complicated in the case of the Breit interaction 
because it is represented by a two-electron operator. General justification of the possibility to simulate the Breit effects 
by means of an one-electron (G)RECP operator with good accuracy and the scheme of such GRECP generation are 
presented in [29] (see also section IC). This scheme is applied in the present work to generate GRECPs for the 
uranium, plutonium, eka-mercury (E112), eka-thallium (E113) and eka-lead (E114) atoms. The 32, 34, 20, 21 and 
22 electrons are explicitly treated in calculations with these GRECPs, correspondingly. Moreover, the 52 electron 
GRECP (52c-GRECP) version for El 12 was also generated. The conventional Coulomb operator for two-electron 
interactions and the point nuclear model should be used in these GRECP calculations. However, they will account 
for the Fermi nuclear charge model that is close to the experimental distribution. Moreover, the Breit interactions of 
the electrons from the state used for the GRECP generation with the electrons explicitly treated in the subsequent 
calculations are simulated by the GRECP (in some sense, the Breit interaction is "frozen" here). 



3 



A. Self-consistent GRECP version for d- and /-elements 

The Self- Consistent (SfC) (G)RECP version [19, 23, 24, 27] allows one to minimize errors for energies of transitions 
with the change of the occupation numbers for the OuterMost Core (OMC) shells without extension of space of 
explicitly treated electrons. It allows one to take account of relaxation of those core shells, which are explicitly 
excluded from the GRECP calculations, thus going beyond the frozen core approximation. This method is most 
optimal for studying compounds of transition metals, lanthanides, and actinides. Features of constructing the self- 
consistent GRECP are: 

1. The all-electron HFDB calculations of two generator states with different occupation numbers N\ and N 2 of the 
OMC d or / shell are carried out for an d- or /-element. 

2. The GRECP versions with separable correction U* 1 and U* are constructed for these generator states em- 
ploying the standard scheme [17, 18, 19, 24]. The GRECP operator with the separable correction has the 
form 



where 



L ' + 1/2 
1=0 j=\l-l/2\ 



- V 



n r , Ij j ' 



(1) 



p y = E l^' m )(^' m l' 



P £«= E |K/im)^)(Kym)^|, 

m=-j 

\lj m ){lj m \ is the projector on the two-component spin-angular function xzj m , |(n c /jm) Ari )((n c ^'m) Ari | is the 
projector on the outer core pseudospinor PnlijXijm, ^nlij an d ^njj are * ne ra dial components of the GRECP 
derived for valence tp^jj an d outer core p^jj P se udospinors for the OMC d or / shell occupation number 
Ni (i=l,2), E^ rc is the core energy, L is one more than the highest orbital angular momentum of the inner 
core spinors and J = L + 1/2. The separable terms (the second and third lines in Eq. (1)) are added to the 
conventional radially-local RECP operator. These terms take into account the difference between the potentials 
acting on the outer core and valence electrons with the same I and j. 

3. The self-consistent GRECP, U sfc , with the quadratic correction writes as 

Sfc _ U^+U^ U N i-U N * f AT N 1 +N 2 \ 
2 + N,-N 2 V 2 ) 



+ B^N omc -^l) , (2) 

where B is some adjustable parameter, iV omc =(^|N omc |'I'), ^ is the many-electron wavefunction for the calcu- 
lated state, and N omc is the occupation number operator of the considered d (/) shell that is written as 

1+1/2 j 

Nome = E/ E/ ^omciomcJin^^t'onJ" 1 ' ^) 

j=\l-l/2\ m=-j 
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i jm (^"omciomcjm) i s the creation (annihilation) operator for the electron in the pseudostate \n omc l omc jm) 
corresponding the original one-electron state \n omc l omc jm), n omc and l omc are the principal and orbital quantum 
numbers of the OMC shell. 

4. The ~Pnjj projectors in U^' from Eq. (2) are replaced by the projectors 

P£u = E \(njjmr)((njjmr\ 

m=-j 

for simplicity, where | (n c /jm) av )((n c Zjm) av | is the projector on the outer core pseudospinor ^jjXijm, 



$Zij (r) = Gno rm [^y (r) + (r)] , (4) 

and Cnorm is the normalizing factor. 

The comparison of self-consistent and conventional GRECP versions by accuracy in calculations on the uranium and 
plutonium atoms can be found in paper [29]. 



B. Term-splitting correction for d, /-elements 



The self-consistent (G)RECP correction gives no improvement in description of splittings to terms, e.g., of the 
configuration 5/| /2^3/2^ s i / 2 of uranium as compared to the parent (G)RECPs [19, 24]. Analysis of the corresponding 
errors shows that the main contribution (about 90 %) is due to smoothing the original OMC spinors in the core 
region. The simplest way to minimize these errors is to use such (G)RECPs, in which the 5/ shell is described by 
nodal pseudospinors, whereas the 4/ pseudospinors are nodeless. To reduce computational efforts, the 4/ shell can 
be treated as "frozen" using the level-shift technique [19, 31]. 

If the small magnitude of the OMC shell (5/ here) relaxation is taken into account, there is another way out that 
can be optimal for the low-lying states. It was suggested in [19] to add the Term-Splitting (TS) correction (see also 
[24]) to the (G)RECP operator 

UTS = ^2 X XlX2 , X3Xi \x 1 )\x 3 )(x 2 \(x4\ 

Xl,X2,X3,X4 

)S WX3 \x 1 )(x 2 \ , (5) 

W Xl,X2,X3 

where X Xl x 2 ,x 3 x 4 is the difference between the two-electron integrals calculated with original spinors and pseudospinors 
for the generator state, the indices w = (n occ l occ j occ m occ ) correspond to the occupied spinors for the calculated state, 
the indices x = (n omc l omc jm) run over all possible j = \l mc ± 1 /2 1 and m = —j, —j + 1, . . . j for the given OMC 
shell. These terms correct the one- and two-electron integrals containing only the 5/ pseudospinors of uranium in the 
considered case. 



C. Accounting for the Breit interaction between different shells 

Let us analyze contributions of the Breit interaction between electrons from different shells to the energy of a heavy 
atom [27]. We will use the estimate (e.g., see [32]) 

{P,P'\{a v a v )\P,P') ^ ^{{vp-v pi )) ■ 

for an uncoupled one-electron state P: (P\a\P)^^r- , ^^ ~aZp, where &i are 4x4 Dirac matrices for the ?-th 
electron, c and v arc velocities of light and electron, a~j^ is the fine structure constant. In the above expression a 
"pseudocharge" , Z* P , is introduced which can be most naturally defined in our consideration as [29] 



Z* P = (P\-\P) , 



(6) 
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that coincides with the nuclear charges only for nonrelativistic electrons occupying the ground states in hydrogen-like 
ions. Besides, (^} can be estimated as (i) for the outermost of the one-electron states P,P' [27]: 



(P,P'\ — |P,P'} ~ min 

ri2 



r r 



= min [Z* P ,Z* pl ] 



As a result, the Breit interaction between the one-electron states P and P' can be estimated as 

B PP , « a 2 ZpZ* P , ■ min [Z* P , Z* P ,] ■ T 

where the correcting factor T ~ [0.1 -j- 1] is introduced, which depends on Al=\l P — lpi\, Aj=\j P — jp>\, etc. 

Applying Eq. (6) for inner core (P=/), outer core (P=c) and valence (Fed) electrons one has Z*^3, 
Z*r^l by the order of magnitude (Zp differs from an "effective charge" of the core with respect to the electron in the 
P-th state, Zp f =Z— , that is usually used in RECP calculations, where Z is the nuclear charge, is the number 
of core electrons with respect to the P-th state). Therefore, Bppt=T~ Y Bppi is as 

Bff> ~ 10 000 000 cm" 1 , B fc ~ 9 000 cm' 1 , B fv ~ 1000 cur 1 , 
B cf ~ 9 000 cm' 1 , B cc > ~ 270 cm" 1 , B cv ~ 30 cur 1 , 

B vf ~ 1000 cm' 1 , B vc ~ 30 cm" 1 , B vv > ~ 10 cm" 1 . 

Let us consider approximations in accounting for the Breit interaction, that we made when outer core and valence 
electrons are included in GRECP calculations with Coulomb two-electron interactions, but inner core electrons are 
absorbed into the GRECP. When both electrons belong to the inner core shells, the Breit effect is of the same order 
as the Coulomb interaction between them. Though Bff does not contribute to "differential" (valence) properties 
directly, it can lead to essential relaxation of both core and valence shells. This relaxation is taken into account when 
the Breit interaction is treated by self-consistent way in the framework of the HFDB method [33, 34] . 

The inner core electrons occupy closed shells. The only exchange part of the two-electron Breit interaction between 
the valence, outer core and inner core electrons, Bf v and P/ c , gives non-zero contribution. The contributions from 
Bf v and Bf c , are quite essential for calculation at the level of "chemical accuracy" (about 1 kcal/mol or 350 cm -1 for 
transition energies). This accuracy level is, in general, determined by the possibilities of modern correlation methods 
and computers already for compounds of light elements. Note, that the contribution from the exchange interaction 
is not smaller than that from the Coulomb part [29]. The inner core electrons can be considered as "frozen" in most 
physical-chemical processes of interest. Therefore, the effective operators for Bf v and Bf c acting on the valence and 
outercore shells, Bj* and Bj*, are of the same kind as the exchange f—v and f—c contributions of the SCF field in 
the Huzinaga-type potential, i.e. these terms can be well approximated by the spin-dependent potential of the form: 



B fv + B % = E ^fMPu + E W?M r ) - ^f»] Pn e «, 

Ij njj 

which has basically the same spin-angular structure as the GRECP has. Thus, it can be taken into account directly 
when the HFDB (not HFD) calculation [8] is performed to generate outer core and valence bispinors but in the 
inversion procedure of the HF equations for generating the components of GRECP, the conventional interelectronic 
Coulomb interaction should be used instead of the Coulomb-Breit one. Then, in the GRECP calculations one should 
consider only the Coulomb interaction between the explicitly treated electrons. 

Due to small relaxation of outer core shells in most processes of interest, these shells can be also considered as 
"frozen" when analyzing the Breit contributions and the P cc < and B cv terms can be taken into account similarly to 
the Bf c and Bf v ones. The error of this approximation will be additionally suppressed by relative weakness of the 
Breit interaction with the outer core electrons as compared to the inner core ones. We note here, that the estimates 
for Z*, Z* and, therefore, for B cc > 7 B cv and B vv > given above are rather the upper limits. For heavy atoms these 
Breit contributions are smaller approximately by one-two orders of magnitude. This decrease is due to enlarged radii 
of the valence and outer core shells and other effects in heavy atoms [29]. For example, for uranium (Z = 92) one has 
Zj-Jnonrel. SCF]-92.4, Z 1 * S [DHFB]~122.4 (starting from Z ~ 30, Z\ s grows faster than Z due to relativistic effects, 
whereas is essentially smaller than the corresponding effective charge Z^[ for all other nl), Z^^l, Zg s ~l, Z£ p ~0.7, 
Zg d ~0.4, Zy s ~0.3. Thus, B cc ', B cv , and B vv i contributions are negligible for the "chemical accuracy" of calculation. 
Therefore, the above made estimates provide us a good background for approximating the Breit interaction by a 
one-electron GRECP operator that should work well both for actinides and for superheavy elements. The numerical 
tests of the GRECPs accounting for the Breit effects are discussed in the next section. 
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II. RESULTS AND DISCUSSION 

For all-electron calculations, we used the atomic HFDB code [35, 36] which allows one to account for the Brcit 
interactions both in the framework of the first-order perturbation theory (PT-1) and by the self-consistent way as well 
as to account for different models of nuclear charge distribution. For test calculations with (G)RECPs, the atomic 
Hartree-Fock code in the jj-coupling scheme (hfj) [17] was used (that was quite sufficient for studying errors of the 
one-electron (G)RECP operators). Both the codes are numerical that allows us to exclude the errors due to the 
incompleteness of basis sets when estimating accuracy of different RECPs and GRECPs. 

The transition energies between states averaged over the low-lying configurations of SHEs 112, 113, 114 and actinides 
U, Pu are presented in tables I, II, III and IV, V, respectively. One can see that the errors due to the point nuclear 
model reach a few thousand wave numbers for the SHEs and several hundred wave numbers for the actinides. The 
considered small variations in the nuclear charge distribution (including the nuclear size) in the framework of finite-size 
nuclei lead to change of the transition energies for the studied SHEs less than on 60 cm -1 . The differences between 
the results with the PT-1 and self-consistent ways of accounting for the Breit interaction are within 7 cm -1 for SHEs 
and actinides whereas neglecting the Brcit effects leads to the errors up to a few thousand wave numbers for the 
studied actinides and several hundred wave numbers for the SHEs. 

The GRECP errors in reproducing the results of the all-electron HFDB calculations with the Fermi nuclear model 
are collected into two groups. First, the GRECP errors for transitions without change in the occupation number of 
the 6c? shell for the SHEs (tables I, II and III) and the 5/ shell for the actinides (tables IV and V) are relatively small 
whereas the corresponding errors of the other tested RECPs for the SHEs are significantly higher. The same number 
of electrons is explicitly treated in calculations with different (G)RECP versions for a given atom. Here and further, 
we do not discuss the particular case of the 52e-GRECP for El 12 if the opposite is not explicitly stated. The RECPs 
of other groups for uranium were tested in paper [19]. It should be noted that they do not take into account the 
large contribution from the Brcit interaction. The Breit effects were also not considered at the generation stage of the 
RECP of Nash et al. [37]. However, it can not explain the large errors for this RECP in tables I, II and III. It is not 
clear from paper [37] which nuclear model was used there. The Breit interaction was taken into account only in the 
PT-1 approximation at the generation stage of the PseudoPotential (PP) of Seth et al. However, the corresponding 
changes in the transition energies are negligible in comparison with the PP errors. 

Second, the GRECP errors for transitions with excitation of one 6d electron for the SHEs or one 5/ electron for the 
actinides are about 400 cm" 1 . These errors have a systematic nature (unlike the corresponding errors for the tested 
RECPs of other groups) and are connected with the fact that the OMC 6d shell for the SHEs and the OMC 5/ shell 
for the actinides in the present GRECP versions are described with the help of nodeless pseudospinors. Obviously, 
these errors can be reduced significantly if one includes the 5d, 5/ electrons for the SHEs and the 4/ electrons for the 
actinides explicitly in the GRECP calculations (see the 52c-GRECP results for E112 in table I). The corresponding 
pseudospinors can be then "frozen" in these GRECP calculations with the help of the level-shift technique [19, 31] to 
reduce the computational efforts. Alternatively, the self-consistent GRECP method described in section I A can be 
used. 

The energies of splittings between terms are considered in table VI for El 12 and table VII for U. The errors of the 
RECP and GRECP approximations and the errors caused by neglecting the Breit effects are within 200 cm" 1 for E112 
(except for the RECP of Nash et al.). The Brcit contributions to the term-splitting energies for U are within 100 cm" 1 
whereas the GRECP errors are up to 750 cm -1 . The latter can be reduced drastically by applying the term-splitting 
correction (see section IB and table VII). The results show that addition of the term-splitting correction allows 
one to reduce the most serious errors up to 10 times for the splittings into terms, thus reducing the errors for the 
energies of transition between terms to the same order of magnitude as the errors for transitions between the states 
averaged over the configurations (when only the self-consistent GRECP is applied) . Obviously, any transition between 
two different terms having different occupation numbers of the OMC shell, N 1 mc and N^ mc , can be presented as a 
combination of three consequent transitions: transiton from the first term to the average over the configuration with 
the same Nj mc , transition between averages over configurations with N5 mc and N^ mc and transition from the latter 
to the second term with N^ mc . Therefore, applying of both the self-consistent and term-spitting GRECP corrections 
to treatment of transitions between any terms allows one to reduce dramatically the (G)RECP approximation errors 
without increasing the number of explicitly treated core electrons of a considered d, /-element. 

In tables VIII and IX, the matrix elements of < r 2 > and radial integrals 2 J R \ fn v ij(r)[f nv ij(r) — tp n jj{r)] \ dr 
(where / n „y is the large component of the Dirac spinor, tp nv lj is the radial part of the corresponding pseudospinor 
and R n is the radius of the last spinor node) are considered for the cases of spinors from different configurations of 
El 12. The errors in these matrix elements and integrals characterize the quality of reproducing the electronic density 
in outer core and valence regions of the atom. One can see that the GRECP allows one to reproduce the electronic 
density in the valence region (the ls\ji and 7p\/2 spinors) with very high accuracy. The one-electron energies for 
spinors from different configurations of El 12 are presented in table X. Similar conclusion can be made in the latter 
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case. 



Conclusions 



Different nuclear models and contributions of the Breit interaction between valence, inner and outer core shells of 
uranium, plutonium and superheavy elements E112, E113, and E114 are considered in the framework of all-electron 
four-component and (G)RECP methods. It is concluded on the basis of the performed calculations and theoretical 
analysis that the Breit contributions with inner core shells must be taken into account in calculations of actinide and 
SHE compounds with "chemical accuracy" whereas those between valence and outer core shells can be omitted. 

The differences in the atomic energies between the cases of the PT-1 and self-consistent ways of treating the Breit 
interaction as well as small variations in the nuclear charge distribution in the framework of finite-size nuclei are not 
essential for the considered accuracy of calculations. However, the difference between the point and finite nuclear 
models is important for the valence (transition) energies. The effects of accounting for the Breit interaction and 
finite nuclear model can be simulated by GRECPs with very good accuracy when only Coulomb interaction between 
the explicitly treated electrons is taken into account. Thus, the GRECP method allows one to carry out reliable 
calculations of actinides, SHEs and their compounds at the level of "chemical accuracy" . 
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TABLE I: Transition energies (TE) between states averaged over the relativistic configurations of E112 derived from HFDB 
calculations with Fermi nuclear model and the corresponding absolute errors of all-electron and (G)RECP calculations (in cm -1 ). 





HFDB HFDB 


HFDB 


HFDB HFD+B 


HFD 


52e- 


20e- 


Ionic 20e-RECP 20e-PP 




(Fermi, 


(Ball, 


(Fermi, 


(Point) 


(Fermi, 


(Fermi, 


GRECP GRECP 


20e- 


of Nash of Seth 




A=296) A: 


=296) A=285) 




A=296) A=296) 






RECP 


et al. 


et al. 




(a) 


(b) 


(a) 


(c) 


(d) 


(e) 


(f) 


(f) 


(g) 


00 


0) 


Configuration 


TE 










Absolute errors 










6dy 2 6d e 5/2 7s 2 1/2 -> 
























6rf3 /2 6df /2 7s 1/2 7p 1/2 


46406 


-3 


22 


1768 


1 


-27 


1 


-17 


588 


3198 


153 


6^3/2 6^5/2 7*1/2 7^3/2 


64559 


-4 


25 


1964 


-1 


239 


4 


-29 


820 


5480 


27 


6^3/26^5/2 7s 1 / 2 8s 1 / 2 


72571 


-3 


22 


1760 


-1 


257 


6 


-25 


719 


5085 


105 


64 /2 6d« /2 7 S i /2 74 /2 


81845 


-4 


23 


1879 


-1 


277 


6 


-18 


809 


5465 


99 


6d% /2 6dl /2 7sl /2 7p{ /2 


28701 


1 


-8 


-644 


2 


-576 


31 


305 


-422 


-3723 


380 


&dl /2 Qdl /2 7s\ /2 7p} i/2 


52595 


1 


-6 


-464 





-267 


37 


277 


-181 


-1254 


189 


&d\/ 2 6d\i 2 7s\/ 2 &s\i 2 


62635 


2 


-10 


-776 





-252 


43 


314 


-315 


-1879 


326 


6d 3/ , 2 6d 5 / 2 7s 1 / 2 7d 3 ^ 2 


72443 


1 


-9 


-666 





-234 


43 


322 


-220 


-1514 


314 


6d3/ 2 6c?g/27si/2 


84449 


1 


-9 


-672 





-234 


43 


322 


-224 


-1531 


308 


6^3/2 6^5/2 7s?/ 2 7Pi/ 2 


53581 


2 


-10 


-765 


2 


-281 


45 


387 


-376 


-3903 


22 


6^3/2 6d 5 / 2 7s 1 / 2 7p 3 / 2 


75273 


1 


-8 


-600 





7 


52 


437 


-84 


-1515 


-126 


6g! 3 / 2 6g? 5 /2 7s 1 / 2 8s 1 / 2 


85677 


2 


-12 


-915 


-1 


25 


60 


477 


-213 


-2126 


22 


6d a 1 2 6d 5 / 2 7s ! / 2 7d 3 / 2 


95546 


2 


-10 


-805 


-1 


43 


60 


484 


-119 


-1760 


7 



(a,c) All-electron Hartree-Fock-Dirac-Breit (HFDB) calculations with Fermi and point nuclear charge distributions, 
accordingly. 

(b) All-electron HFDB calculation with the uniform nuclear charge distribution within a sphere. 

(d) All-electron HFD calculation with accounting for the Breit interaction within PT-1 (HFD+B) and with Fermi nuclear 
model. 

(e) All-electron HFD calculation without accounting for the Breit interaction (HFD) and with Fermi nuclear model. 

(f) GRECP generated in the present work from HFDB calculation with Fermi nuclear model. 

(g) Semi-local RECP generated here from HFDB calculation with Fermi nuclear model on the ionic closed-shell 
generator-state. 

(h) RECP from [37] generated from HFD calculation. 

(i) PP from M. Seth et al. to be published (P. Schwerdtfeger, private communication, 2003) generated from HFD+B 
calculation. 



[30] N. S. Mosyagin and A. V. Titov, arXiv.org/ physics/0406143 (2004); J. Chem. Phys., in press (2005). 

[31] A. V. Titov, N. S. Mosyagin, A. B. Alekseyev, and R. J. Buenker, Int. J. Quantum Chem. 81, 409 (2001). 

[32] L. N. Labzowsky, G. L. Klimchitskaya, and Y. Y. Dmitriev, Relativistic Effects in the Spectra of Atomic Systems (Institute 

of Physics Publishing, Bristol and Philadelphia, 1993), 340 pp. 
[33] H. M. Quiney, I. P. Grant, and S. Wilson, J. Phys. B 20, 1413 (1987). 

[34] E. Lindroth, A.-M. Martensson-Pendrill, A. Ynnerman, and P. Oster, J. Phys. B 22, 2447 (1989). 
[35] V. F. Bratzev, G. B. Deyneka, and I. I. Tupitsyn, Bull. Acad. Sci. USSR, Phys. Ser. 41, 173 (1977). 

[36] I. I. Tupitsyn and A. N. Petrov, in 5-th Session of the V.A. Fock School on Quantum and Computational Chemistry 

(Novgorod the Great, 2002), p. 62. 
[37] C. S. Nash, B. E. Bursten, and W. C. Ermler, J. Chem. Phys. 106, 5133 (1997), [Erratum: JCP 111 (1999) 2347]. 
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TABLE II: Transition Energies (TE) between states averaged over the relativistic configurations of E113 (in cm 1 ). See Table 
I. 





HFDB HFDB 
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HFDB HFD+B 


HFD 




Ionic 21c-RECP 21e-PP 
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A=297) A 


=297) A=284) 
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et al. 


et al. 
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TE 
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6dS /a 6dg /2 7fl? /2 7pJ /a 
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3 
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-2 
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-21 
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-349 


6d3/26d5/ 2 7si /2 8si /2 


34981 





-2 


-128 


-3 


354 


5 


112 


-205 


-307 


60(3/26^5/27*1/27^3/2 


45172 








-4 


-3 


374 


9 
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140 


-275 


6^3/26(^5/2 7si/ 2 6/5/2 


50338 
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-3 


374 


9 


196 


127 


-276 


6^3/2 6c?5/ 2 7Si/ 2 5(77/2 
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-10 


-3 


374 


9 
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127 


-276 
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-3 
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9 
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-276 


6d3/ 2 6dg/ 2 7s}/ 2 7pi/ 2 
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-4 


32 


2220 


2 


-60 


28 
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4830 
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6d3/26d 5 / 2 7si/27pi/27p 3 /2 
6^ /2 6dg /2 7ai /2 7^ /a 


83184 


-5 


36 


2485 


-1 


241 


-6 


833 


5170 


-172 


112678 


-6 


41 


2843 


-3 


612 


-10 


1171 


5717 


-504 


6^3/261^5/2 7si/ 2 7pi/ 2 


115758 


-5 


34 


2344 


-1 


250 


-3 


784 


5143 


-105 


6d3/ 2 6dg/ 2 7si/ 2 7p3/2 


149550 


-5 


40 


2739 


-3 


654 


-9 


1163 


5784 


-454 


6d3/ 2 6dg/ 2 7si/ 2 


234435 


-5 


37 


2583 


-4 


747 


-2 


1221 


6102 


-336 


6^3/26^5/2 7si /2 7Pi/2 


47371 


2 


-13 


-864 


3 


-739 


404 


-597 


-2349 


322 


6^3/266(5/2 7Sl/27pi/27p3 /2 


74898 


1 


-9 


-606 


1 


-391 


344 


-378 


-2055 


-44 


6^3/26^5/2 7s? /2 7p3/ 2 


110406 


1 


-4 


-244 


-2 


22 


310 


-47 


-1528 


-407 


6o!3/ 2 6d5/27si/ 2 7pi/ 2 


110120 


2 


-13 


-882 





-388 


386 


-451 


-2298 


41 


6o!3/26d 5 / 2 7si/27p3/2 


150102 


1 


-7 


-477 


-2 


59 


344 


-82 


-1667 


-339 


6^3/260(5/2 7Si/ 2 


239523 


2 


-12 


-807 


-2 


144 


416 


-39 


-1617 


-188 


6d 3 3/2 6d 6 5/2 7s 2 1/2 7p 2 1/2 


78821 


2 


-15 


-983 


2 


-375 


380 


-649 


-2230 


-270 


6dl /2 6dl /2 7s 2 l/2 7p[ /2 7pl /2 


104059 


1 


-11 


-742 





-49 


412 


-364 


-1931 


-544 


6d 3 / 2 6d 6 / 2 7s 1 / 2 7p 3 / 2 
6d 3 3/2 6di /2 7s 2 1/2 7p{ /2 


137048 


1 


-6 


-403 


-2 


341 


481 


38 


-1402 


-804 


139819 


2 


-15 


-1021 





-42 


447 


-439 


-2161 


-466 


6dl /2 6dl /2 7s 2 1/2 7pl /2 


177137 


1 


-9 


-638 


-3 


381 


516 


9 


-1523 


-736 
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TABLE III: Transition Energies (TE) between states averaged over the relativistic configurations of E114 (in cm x ). See Table 
I. 





HFDB HFDB 


HFDB 


HFDB HFD+B 


HFD 




Ionic 22e-RECP 22e-PP 
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(Point) 
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(Fermi, 22e-GRECP 22e-RECP 
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of Seth 




A=298) A 


=298) A=289) 




A=298) A=298) 






et al. 


et al. 


Configuration 


TE 










Absolute errors 










6d 3 ^ 2 6d^ 2 7s 2 ^ 2 7p 2 ^ 2 — > 
&dt /2 &dt /2 7s\ /2 7p\ /2 7p\ /2 






















29093 


-1 


3 


314 


-2 


380 


-46 


211 


449 


-457 


&dt /2 Gdl /2 7s\ /2 7p\ /2 %s\ /2 


41211 





-1 


-135 


-3 
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-51 


53 


-333 


-348 


6d 3/2 6d 6 5/2 7s 2 1/2 7p{ /2 8p{ /2 


48149 








-5 


-3 
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-44 


114 


11 


-318 


Qd\ /2 6dl /2 7s 2 1/2 7p\ /2 7d\ /2 


52230 








16 


-3 
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-44 
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115 


-320 


6dy 2 Qd®/ 2 7s\/ 2 7p\/ 2 Qfl/ 2 


57618 








7 


-3 


384 


-43 


143 


86 


-318 


6d 3 / 2 6d 5 / 2 7s 1 / 2 7p 1 / 2 5g 7 / 2 


60094 








7 


-3 


384 


-43 


143 


86 


-317 


Gd 3 / 2 6d 5 ^ 2 7s 1 / 2 7p 1 / 2 
6d 3 / 2 6d 6 / 2 7s 1 / 2 7p 3 / 2 
6dy 2 6d b 5/2 7s 2 1/2 7pl /2 8sl /2 


64483 








7 


-3 


384 


-43 


143 


86 


-317 


66669 


-2 


6 


755 


-5 


833 


-52 


535 


1209 


-926 


81879 


-1 


3 


277 


-5 


850 


-64 


377 


374 


-834 


6d 3/2 6dt /2 7sl /2 7pl /2 


106776 


-1 


4 


497 


-5 


872 


-53 


503 


969 


-808 


/-> j4 /? j6 h 2 q 2 
Dct 3/2 Dffl 5/2 ' s l/2 8s l/2 


108893 


1 


-3 


-361 


-6 


883 


-82 


174 


-725 


-718 


6d 3/2 6dt /2 7sl /2 8s[ /2 


136567 





-1 


-207 


-6 


907 


-72 


285 


-241 


-680 


bd 3/2 bd 5/2 7s 1/2 


197486 





1 


128 


-6 


yoi 


-45 


547 


853 


-584 


&d\i 2 Qd\/ 2 7s\i 2 7p\i 2 7p\i 2 
&d 4 3/2 &d% 2 7s\ /2 7p\ /2 8s\ /2 


102896 


-6 


24 


3110 





256 


96 


929 


6650 


-327 


115405 


-5 


21 


2745 


-1 


224 


59 


754 


5987 


-243 


Vol 1ir2 

Da3^/ 2 t)a 5 y2 ' s l/2 'Pl/2 


138842 


-6 


23 


2905 


-1 


233 


73 


848 


6439 


-206 


6dt /2 6dl /2 7s 2 1/2 7pl /2 7pl /2 


97736 


2 


-6 


-771 


1 


-506 


472 


-631 


-3156 


-28 


6dt /2 6dl /2 7sl /2 7pf /2 8s[ /2 


112486 


3 


-10 


-1277 


1 


-543 


473 


-830 


-4091 


83 


6d 3/2 6dl /2 7s 2 1/2 7p 2 1/2 


136356 


2 


-9 


-1129 


1 


-534 


487 


-732 


-3647 


122 


Qdl /2 Qdl /2 7s 2 1/2 7p\ /2 7pl /2 
&dl /2 Qdl /2 7s\ /2 7p\ /2 8s\ /2 
Gdl /2 Gdl /2 7s\ /2 7p 2 1/2 


133837 


2 


-7 


-904 





-91 


391 


-746 


-2821 


-675 


149162 
173108 


3 
3 


-11 
-10 


-1415 
-1265 






-126 
-117 


380 
391 


-952 
-855 


-3758 
-3309 


-579 
-541 
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TABLE IV: Transition Energies (TE) between states averaged over the nonrelativistic configurations of uranium (in cm x ). 
See footnotes in Table I. 
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5 f6d 4 


53806 


230 





196 


-15 


5/ 3 7s 2 6rf 1 - 












5/ 5 


99459 


252 


4 


1126 


-671 



12 



TABLE V: Transition Energies (TE) between states averaged over the nonrelativistic configurations of plutonium (in cm x ). 
See footnotes in Table I. 
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TABLE VI: Transition Energies (TE) between terms of E112 (in cm 1 ). See footnotes in Table I. 
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TABLE VII: Transition Energies (TE) between terms of uranium (in cm 1 ). See footnotes in Table I. 
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11 


134 


-29 


139 


-42 


J=5 


2785 





22 


117 


-24 


121 


-36 


J=6 


4606 


1 


29 


181 


-42 


188 


-62 


J=7 


-6030 


1 


78 


-176 


12 


-186 


26 


T O 


-5542 


1 


90 


-149 


6 


-158 


17 


u t /5/2 uti 3/2 ' *l/2 u u 
J=l 


-19109 





23 


-432 


-61 


-426 


-8 


J=2 


-15310 





1 


-304 


-45 


-288 


5 


J=3 


-23656 





41 


-598 


-77 


-593 


-8 


J=4 


-26013 





21 


-638 


-69 


-618 


23 


J=5 


-32544 





36 


-754 


-86 


-732 


21 


J=6 


-39562 





-2 


-724 


-57 


-671 


82 



(a) Term- Splitting (TS) correction generated in the present work from HFDB calculation with Fermi nuclear charge 
distribution. 

(b) Self-Consistent Generalized Relativistic Effective Core Potential (SfC GRECP) generated in [29] from HFDB calculation 
with Fermi nuclear charge distribution. 
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TABLE VIII: Matrix Elements (ME) of < r 2 > for some spinors from states averaged over the relativistic configurations of 
E112 (in a.u.). See footnotes in Table I. 





HFDB 


HFD 


52e- 


20e- Ionic 20e- 


20e-RECP 20e-PP 




(Fermi, 


(Fermi, 


GRECP GRECP 


I V 1 A. I 


of Nash of Seth 




A=296) A=296) 








61 di. 


ei at. 


Config., spinor 


ME 






Absolute 


errors 






64/2 64/2 7^>i/2 
















64/2 


3.150 


n nnK 
-0.005 


n aai 
0.001 


n no /i 
0.024 


0.030 


n r\aa 
O.Udd 


n n70 


64/2 


3. (ai 


n nno 
0.002 


n nm 
U.UU1 


n no a 
U.Uz4 


n nQO 


U.U / 4 


n nccT 


T.Si /<■> 

1 °l/2 


7.157 


-0.023 


n nnn 
0.000 


n nnK 
0.005 


n nnn 

-u.uyy 


n a oc 
-0.425 


n no /i 
0.024 


64/264/2^1/2 


7p{/2 














64/2 


3.144 


-0.004 


0.001 


0.023 


0.028 


0.064 


0.071 


64/2 


3.648 


0.002 


0.001 


0.024 


0.031 


0.069 


0.057 


7S1/2 


6.898 


-0.020 


0.000 


0.002 


-0.097 


-0.394 


0.022 


7P1/2 


13.023 


-0.116 


-0.001 


0.005 


-0.131 


-0.841 


0.055 


64/2 64/2 7«l/2 


7Pl/2 














64/2 


3.057 


-0.004 


0.001 


0.023 


0.030 


0.071 


0.071 


64/2 


3.522 


0.002 


0.001 


0.025 


0.035 


0.080 


0.059 


7S1/2 


6.739 


-0.019 


0.000 


0.001 


-0.092 


-0.361 


0.025 


7P1/2 


11.259 


-0.087 


-0.002 


-0.001 


-0.105 


-0.597 


0.049 



TABLE IX: Radial integrals 2 J~ | f nM (r)[/„„y(r) — V 5 n t ,!i( r )] I ^ r f° r valence spinors from states averaged over the relativistic 
configurations of E112 (in a.u.) where fn v ij is the large component of the Dirac spinor from HFDB calculation with the Fermi 
nuclear charge distribution for A — 296, <fn v ij is the radial part of the corresponding pseudospinor (or the large component of 
the Dirac spinor), R n is the radius of the last node for the spinor. See footnotes in Table I. 





HFD 


52e- 


20e- Ionic 20e- 


20e-RECP 20e-PP 




(Fermi, 


GRECP GRECP 


RECP 


of Nash 


of Seth 




A=296) 








et al. 


et al. 


Configuration, spinor 






Integrals 






6^3/26^5/27^1/2 














7S1/2 


0.0037 


0.0000 


0.0006 


0.0131 


0.0590 


0.0024 


64 4 /26^/ 2 7si/ 2 7pi/ 2 














7«l/2 


0.0036 


0.0001 


0.0002 


0.0137 


0.0581 


0.0022 


7P1/2 


0.0087 


0.0001 


0.0004 


0.0091 


0.0610 


0.0047 


64 4 /264 5 /27s 2 / 2 7pi/ 2 














7S1/2 


0.0034 


0.0001 


0.0003 


0.0134 


0.0546 


0.0026 


7P1/2 


0.0079 


0.0002 


0.0002 


0.0088 


0.0514 


0.0050 
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TABLE X: One-electron energies, e, for some spinors from states averaged over the relativistic configurations of El 12 (in a.u.). 
See footnotes in Table I. 





HFDB 


HFD 


52e- 


20e- Ionic 20e- 


20e-RECP 20e-PP 




(Fermi, 


(Fermi, 


GRECP GRECP 


RECP 


of Nash 


of Seth 




A=296) A=296) 








et al. 


et al. 


Config., spinor 


£ 






Absolutc 


errors 






6^3/26^5/2 7s 1/2 

6^3/2 


0.5624 


0.0003 


0.0000 


0.0001 


-0.0026 


-0.0112 


-0.0037 


6^5/2 


0.4432 


-0.0011 


0.0000 


-0.0001 


-0.0026 


-0.0095 


-0.0009 


7S1/2 


0.4497 


0.0014 


0.0000 


-0.0003 


0.0042 


0.0272 


0.0001 


6rf3/26dg / / 2 7sj;/ 2 

6^3/2 


7rf/2 

0.6148 


-0.0001 


0.0000 


-0.0001 


-0.0025 


-0.0118 


-0.0039 


6^5/2 


0.4870 


-0.0017 


0.0000 


0.0000 


-0.0024 


-0.0108 


-0.0009 


7S1/2 


0.5217 


0.0011 


0.0000 


0.0000 


0.0048 


0.0278 


0.0004 


7P1/2 


0.2248 


0.0015 


0.0000 


-0.0001 


0.0012 


0.0114 


-0.0004 


6c?3 / , 2 6dg / , 2 7sJ //2 

6^3/2 


7Pi/ 2 
0.6663 


-0.0002 


0.0000 


-0.0002 


-0.0042 


-0.0173 


-0.0043 


6ds/2 


0.5314 


-0.0018 


0.0000 


-0.0001 


-0.0040 


-0.0159 


-0.0014 


7«l/2 


0.5253 


0.0010 


0.0000 


0.0000 


0.0044 


0.0258 


0.0001 


7P1/2 


0.2653 


0.0017 


0.0000 


0.0000 


0.0011 


0.0110 


-0.0005 



